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Abstract We will comment on the perihelion/periastron advance of celestial bodies due to 
the cosmological constant A. It is well known that the cosmological constant A causes 
the perihelion/periastron shift; however, there seems to still exist a discrepancy among 
the various derived precession formulae. We will point out that the expression Aoa = 
{nc^Aa^ /{GM))y/\-e^ is the general formula for any orbital eccentricity e and the ex- 
pression A(Oa = {nc-Aa^ / {GM)){\ — e^)^ comes from the nearly circular (e <C 1) approx- 
imation. 

Keywords Celestial Mechanics ■ Gravitation ■ Cosmological Constant/Dark Energy 
1 Introduction 

The 201 1 Nobel Prize in physics was awarded for the discovery of the accelerating expan- 
sion of the Universe |[T]l2l. According to the present literature, theory that is most suitable 
for explaining this phenomenon is that of the cosmological constant A , generally referred 
to as dark energy. Although the cosmological constant or dark energy can be considered to 
account for the current expansion of the Universe in the simple terms, the details are still far 
from clear. Meanwhile, several attempts have also been made for explaining this accelerating 
expansion without dark energy ll3l l4ll5ll6ll7ll8ll9lfT0l . 

To verify the existence of the cosmological constant/dark energy and/or the potential- 
ity of alternative gravitational theories, it is most natural to investigate their properties in 
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the context of classical tests of general relativity, i.e., the perihelion/periastron advance of 
celestial bodies and the deflection of light. To this end, we must elucidate the difference 
in such effects among the cosmological constant and the alternative gravitational theories. 
Therefore, it is important to derive a rigorous formula for such effects resulting from the 
cosmological constant A . Recently, the effect of the cosmological constant on the bending 
of light rays was intensively investigated and discussed by many authors; the detailed reports 
can be found in llTTnMMM5]fTM T7l [T8]fT9]|20l and the references therein. 

The influence of the cosmological constant on the orbital motion is a historical problem 
(see section 45 of 1211 ). In principle, it is clear that the cosmological constant contributes to 
the perihelion/periastrion advance, although this contribution is too small to detect because 
of A w 10^^^ m^^. However, there is a discrepancy in the obtained precession formulae, in 
terms of the difference in the eccentricity dependence ( 1 — e^)^ l21ll22lfT8l and \/\— 1231 
l24l . The contribution of cosmological constant on the orbital dynamics is also investigated 
in other papers Il25ll26ll27ll28ll29ll30ll31ll32ll33ll34l , and constraints on the cosmological con- 
stant from the Sun's motion through the Milky Way is recently discussed [351. 

In fact, this issue has been resolved in 1361 using the general framework of radial/central 
force perturbation. Subsequently, it was shown that the physical interpretation of the formula 
derived in 1361 is the precession of Hamilton's vector 11371 . A similar problem as that in B6I 
was treated in terms of the Gaussian planetary equation II38I . and Multiple Scales Method 
1391 . In I36II38I[39]| , it was shown that the correct eccentricity dependence is as 
previously shown in 123 11241 . However, in spite of these facts, some confusion still seems 
to exist (Refer to e.g., II18I ). Therefore, the purpose of this short note is to arrange the 
discussions on the perihelion/periastron advance due to the cosmological constant A . 

2 Perihelion/periastron advance due to cosmological constant 

2. 1 Analytical approach by direct integration of perturbation potential/force 

In order to examine the additional perihelion/periastron advance due to A, let us assume the 
spacetime to be spherically symmetric; without loss of generality, we write the metric as 

^/^Z = - [l - ^ + c^dt^ + [l - 7 + W'(^)] + r^dn^, (1) 

where rg = 2GM/c^, dQ~ = dQ- + sin" 

is an additional term that is a func- 
tion of r and generally depends on the gravitational theories. Then, the geodesic equation of 
the test particle with mass m becomes 
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^^-^(l + ^jtl-^+^W]. (3) 

Here, T is the proper time, and E and L are two constants of motion: the total energy and 
angular momentum, respectively. According to the standard approach, we define the New- 
tonian potential 4>(/-) as 



0ir) ^ 1 jirn^ [U{r)-c^] = - ^ + ^ + ^-W{r)c\ (4) 



where we set .if = L/m and let S" = {\/2){E /mc) , then, we have. 



1 f dr^^ 



2UJ +'^>W = ^. (5) 

If we rewrite V{r) = {\ /2)W{r)c~ and T w ? (f is coordinate time), Eq. l|5j becomes com- 
pletely equivalent to Eq. (15) of I36II . Thus, henceforth, we follow the approach employed 
in (361. 

The perihelion/periastron shift Adp due to the perturbation potential V{r) with respect 
to Newtonian potential —GM/r is expressed as (see Eq. (30) in 1361 ) 



_ z& dV{z) 

~ GMe^ J -I yrr^ dz ' 



AG, = / (6) 



in which r = ^ / {I + ez) and £ = .if^/(GM) = a{l — e^), a being the semi-major axis of the 
orbit. It should be noted that in 1361 , V{r), or equivalently V(z), is defined including the 
mass of the orbiting particle m (refer to Eq. (15) in II36I ). However, in our case, V{r) or V(z) 
does not contain m; hence, hereafter, in does not appear in the precession formula. 

We are interested in the power-law perturbing potentials of the type V{r) = Unr". Noting 
the relation, 

^('•) = -«-(„+i)'--f"+'' =-«-(«+i)(l+«)''+V^"+', (7) 
we have the precession formula, 

4,(-(„ + l)) = -—^J_^^_dz. (8) 

For discussing the effect of cosmological constant A, the Schwarzschild-de Sitter or Kottler 
metric l40l . 



ds' = -\^l-^--jrjc'dt'+\^l-^--jr'j dr' + r'dQ' (9) 

is widely used and rg/r term causes famous general relativistic perihelion/periastron shift. 
In the case of perihelion/periastron advance due to A, we consider W{r) = — (l/3)Ar^ and 
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then V{r) = —{l/6)Ac^r^. It should be noted in this case, — («+ 1) = 2; thus, n = —3 and 
«2 = -(l/6)Ac2. 

The integral part in Eq. ([8} is 



1 z{l+ez)-\ -ixey/X-e^ 

Therefore, from Eqs. (|8) and (IIP) , the perihelion/periastron advance due to A is given by 
the following equation (again, note £ = a(l — e^) and n = —3) 



KcrAa 



A(Oa= VI -e^. (11) 

Eq. jl lb is in agreement with the results of Il23ll24ll38l . 

Let us summarize the outline of above derivation. We start from the SdS/Kottler metric 
Eq. ^ and its "Newtonian approximation (weak field approximation)". Then the pertur- 
bation potential Ua = — (l/6)Ar2 or perturbation force Fa = (l/3)Ar is determined. The 
result, Eq. dllt . is obtained exactly by "the direct integration of given perturbation potential 
or force" e.g. Eq. dlOt . with the help of the knowledge on Gaussian hypergeometric function 
2-Fi(o!,jS,7; x), see Eq. (37) in (36). We mention that Eq. dlOb can be shortly calculated by 
the variable transformation. 



COS(^ ■ 



!COS(§ 



(12) 



which produces the integral 



Vl -e2 /■'T 2,71 eVl-e^ 

j^^l icos^-e){l-ecos^)d^ = --j^—^. (13) 

It is worthy to note when integrating Eq. dlOt . the integral part is not expanded with respect 
to the orbital eccentricity e. 



2.2 Standard perturbation method 

The perihelion/periastron advance due to A is considered by means of the standard pertur- 
bation method. Here let us discuss this approach according to 1221 . 

Adopting SdS/Kottler metric Eq. we begin with the second-order geodesic equation 
for time-like world-line, 

d^u GM 3GM 2 1 

^ + "-l|- = ^" -3Zp' " = 7' ^^^^ 
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here L- is the z component of orbital angular momentum^, and right-hand side of l ll4b can be 
considered as the perturbation to the Keplerian motion. The Keplerian motion is described 
by 

and its solution is 

- = M = ^^-(1 -|-ecos(i)). (16) 
r Lt 

Since it is well-known that the first term in right-hand side of Eq. J 14b causes famous general 
relativistic precession formula, Aw^^ = 6TtGM /[ac^{\ — e^)), and Eq. MA\ is linear differ- 
ential equation, then we concentrate on the perturbation due to the cosmological constant 
A; 

(fiu GM Ac^ 

Inserting Eq. ( 116b into right-hand side of Eq. ( 117b . we have 

Ac^ Ac-Lt 1 

(18) 



3Z.?m3 3(GM)3 (l-Fecos(^)3' 

Here, it is worthy to mention when obtaining Eq. (II lb . the part 1/(1 -|-ecos(^)^ is exactly 
integrated with changing the variable cos ^ with z (see Eq. dlOb). However, in the standard 
literature such as ||22j, this part is expanded with respect to the orbital eccentricity e, 

^ l-3ecos(^+6e^cos-(^-Fi^(e^), (19) 



(1 + ecos(^)-^ 

and only ecos(^ term is retained in II22I . Comparing with the derivation process of general 
relativistic precession formula, the precession formula due to A is derived as (seeEq. (14.25) 
in I22ll ) 

nc'-AL*'. 

4co.«^. (20) 

It should be emphasized here although the expansion itself is "truncated at the first order in 
e\ the orbital angular momentum L- is given by 



L- = VGM£ = J GMa{ 1 - e2) (21) 



then the quadratic terms in orbital eccentricity e is remained in the final residt due to L- but 
these terms are not cared in the lowest order perturbation method considered. 



In 1221 . the j component of angular momentum is expressed by h instead of L-. 
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Inserting Eq. (I21t into Eq. ilOl . it is found 



Ttc^Aa^ 2\3 



(22) 



which is also same results of II21II18I . 

If we expand the integral part in Eq. i llOt up to the second order in e and integrate, then 
we have 



1 z{l+ez 



1-3 



/■I z(l-3ez + 6eV) , 3j: 



(23) 



and combining with Eq. ([8]l, Eq. i l22t is recovered. 

The eccentricity dependence of Eqs. i ll lb and (122b is plotted in Fig.[T] It is clear that 
Eqs. dub and ( 122b behave quite differently as the orbital eccentricity becomes large. 



e 



3^ 




0.4 0.6 
Orbital Eccentricity 



Fig. 1 The eccentricity dependence of Eqs. Jilt and i22i . In this plot, we set A « 10 m ^, M = 2.0 ; 
lO^^kg, a= 1.5 X 10" m. 



3 Summary 

In this short note, we commented on the perihelion/periastron advance of celestial bodies 
due to the cosmological constant A. As we stated before, our prescription is based on 11361 , 
in which the problem of perihelion/periastron shift due to the cosmological constant A has 
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practically been resolved. However, we hope that this note contributes to clear up the con- 
fusion about the perihelion/periastron advance due to the cosmological constant A , and is 
helpful in discussing the probability of alternative theories of gravitation since A -like terms 
arise also in various theoretical contexts such as f{T) |41| and /(/?) gravity 1421 . 

Just before closing this letter, we would mention the following issue; for high orbital ec- 
centricity e and large central mass M, the general relativistic value of perihelion/periastron 
advance becomes large, while Eqs. ([77]) and i22i tend to zero (see again Fig.[T]l. This prop- 
erty in Eqs. dllt and l l22t seems to be physically counter-intuitive. As a possibility, this 
counter-intuitive feature may have roots in the fact that essentially, Schwarzschild-de Sit- 
ter/Kottler metric does not become asymptotically flat spacetime. However the derivation 
of Eqs. dlll l and ( 122b may assume implicitly that the background metric is Minkowskian. 
Therefore, the problem on perihelion/periastron advance due to cosmological constant A 
should be investigated and discussed further in a careful manner in the context of real ex- 
panding universe (see, for example, II43I . and similar topic was recently discussed by 1441 '). 
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